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Isoscalar off-shell effects in threshold pion production from pd collisions 
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We test the presence of pion-nucleon isoscalar off-shell effects in the pd — > n'^t reaction around the 
threshold region. We find that these effects significantly modify the production cross section and 
that they may provide the missing strength needed to reproduce the data at threshold. 
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I. INTRODUCTION 

Pion production from nucleon-induced reactions has 
the potential to probe the nuclear phenomena at short 
distance since it involves processes transferring large mo- 
menta to the target nucleus. But the pion also medi- 
ates the nuclear force; hence meson production (or ab- 
sorption) plays a fundamental role in hadron dynamics 
because it may reveal facets of meson-baryon couplings, 
and of meson-exchange processes in general, which would 
remain hidden otherwise. 

In absence of reliable calculations on meson dynam- 
ics within an interacting multinucleon context, one has 
to rely on the determination of the reaction mechanisms 
which dominate the process. Even so, the analysis of the 
process is complicated by the fact that a general treat- 
ment of the reaction mechanisms reveals the occurrence 
of many terms, and one is forced to introduce further 
assumptions in order to reduce the number of terms to 
a few, tractable ones. This reduction clearly introduces 
ambiguities making it more difficult to extract informa- 
tions about the nuclear wave function at short distances, 
or about the modifications of the hadron interactions be- 
cause of the presence of other nucleons. The situation is 
somewhat simplified if we consider nucleon-induced pro- 
duction close to the pion threshold, since there the s-wave 
mechanisms of the elementary NN—^ ttNN inelasticities 
dominate, while the p-wave mechanisms (including the 
isobar degrees of freedom) can be treated as corrections. 

In the past decade, with the aim to clarify the na- 
ture of the elementary NN^ ttNN s-wave inelastici- 
ties, a great deal of experimental and theoretical activity 
has been made in pion production from NN collision at 
energies close to threshold. The situation has been re- 
cently reviewed by Meyer Advances in experimental 
techniques allowed to measure in particular the reaction 
pp TT°pp cross section very close to threshold. This 
reaction filters the s-wave nN coupling in the isoscalar 
channel. Standard theory including the one-body term 
and isoscalar rescattering constrained by the ttN scatter- 
ing lengths underestimated the cross section by a factor 
of five. Unexpectedly, there have been two theoretical 
explanations for this discrepancy, not just one. The en- 



hancement in the cross section can be explained by in- 
voking short-range nucleon-nucleon effects ||^ , where the 
important effects can be simulated by a and uj exchanges 
The explanation is appealing, since in this case the 
pion field couples with the two-nucleon axial charge op- 
erators, and this coupling provides an explicit link to 
the inner part of the nucleon-nucleon interaction, which 
is notoriously difficult to disclose. But the results have 
been entirely explained by resorting also to an off-shell 
enhancement of the isoscalar nN amplitude 

The calculation by Hernandez and Oset employs an s- 
wave rescattering diagram where the isoscalar amplitude 
is described by the combined effect of a strong short- 
range repulsion and a medium range attraction of similar 
strenght, where the repulsion is represented by a contact 
term and the attractive part is parameterized by means 
of a tr-exchange diagram. Originally, this representation 
has been derived with dispersion theoretic methods by 
analyzing the experimental data on pion-nucleon scatter- 
ing with discrepancy functions Q. Subsequently, the fit 
in Ref. |^ has been reinterpreted as being generated by a 
sigma-exchange term plus a w-channel term including N 
exchange (or virtual NN creation) and other short-range 
contributions 

That s-wave pion production/absorption is governed 
by off-shell effects has been known for quite a few years 
1^,^. Hachenberg and Pirner used a field theoretical de- 
scription for ttN scattering in s wave based on the linear 
a model and on pseudoscalar t:NN interaction. This 
combination results in a large cancellation between the 
cr-exchange diagram and the nucleon Born terms. The 
cancellation however breaks down producing an enhance- 
ment when one pion leg is off-mass shell, as happens 
when the pion rescatters before being absorbed. On the 
contrary, in the isovcctor channel the cancellation occurs 
more efficiently off-shell, producing a suppression of the 
charge-exchange interaction. Thus, the relative impor- 
tance of the isospin odd and even channels is reversed 
because of the off-shell extrapolation. Yet, another off- 
shell extrapolation has been considered in Ref. [^, de- 
rived from a field theoretic model of the ttN interaction 
which has been developed at Jiilich This approach 
is similar in spirit to the model developed by Hachenberg 
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and Pirner, but here the meson exchanges in the scalar 
and vector channels are derived from correlated 2tt ex- 
changes. The rescattering diagram of Ref. |^ has a less 
pronounced off-shell enhancement with respect to both 
■nN models of Refs. |10 and and one must introduce 
here two-nucleon short-range mechanisms such as those 
mediated by a and uj exchange in order to reproduce the 
total cross section for pp ppiT° at threshold. Finally, 
the problem of the off-shell extrapolation for the pion 
rescattering mechanism has been considered also in the 
more systematic framework of chiral-perturbation theory 
(xPT) |jri[]. The main feature here is that the rescatter- 
ing diagram, once extrapolated off-shell, produce a neg- 
ative interference with the one-body term, thus yielding 
a cross section substantially smaller than the measured 
ones. In this case the inclusion of heavy meson-exchange 
effect does not solve the discrepancy. However, in another 
chiral-perturbation approach based on full momentum- 
space treatment , the rescattering diagram was shown 
to be larger by a factor of 3, thus leading to a substan- 
tial increase at threshold. More recently, the xPT ap- 
proach has been carried further on by considering pion 
loop diagrams which might simulate cr-exchange phenom- 
ena, with the finding that these higher order contribu- 
tions provide important improvements, but questioning 
at the same time the convergence properties of the power 
counting expansion for the specific reaction under consid- 
eration [ p^ . 

The three-nucleon system is a richer testing ground 
for studies of pion production and scattering. The ad- 
dition of just one nucleon increases the complexity of 
the reaction which involves now the simplest nontrivial 
multinucleon system where it is possible to test the fun- 
damental NN^ ttNN process and, at the same time, 
to solve the accurately the nuclear dynamics with Fad- 
deev methods. Applications of Faddeev methods to pion 
production/absorption started very recently with stud- 
ies centered around the A resonance |lj,|l5| and her cm 
we apply the same technique of Ref. [15| to study pion 
production at threshold. Besides the obvious difhculty 
of performing calculations with three nuclcons instead of 
two, one encounters here the additional difficulty that 
for the pd n'^t reaction one must include from the 
start the effect of p-wave mechanisms, on top of the s- 
wave ones. This contrasts with the findings for the two- 
nucleon case, where the effect of the p-wave mechanisms 
(including the A), tends gradually to zero in approach- 
ing the threshold limit. The effect of this difference can 
be immediately perceived in the differential production 
cross section, since for NN collision the angular depen- 
dence evolves gradually with energy, while in the case 
of pd collisions it exhibits a remarkable s- and p-wave 
interference in the threshold region, with strong forward- 
backward asymmetry ]T6[ |. 

In this work, we have centered our study on the ef- 
fects due to the s-wave rescattering processes, taking 



into account both isoscalar and isovector components and 
their interference to the p-wave mechanism (containing 
also the A degrees of freedom). We have in particular 
taken into account the off-shell effects in the isoscalar 
channel by following the same prescription suggested in 
Ref. [Q to explain the size of the excitation function for 
the pp — + ppT:° process. It is important to stress that 
there are still large uncertainties inherent to this off-shell 
extrapolation, and such calculations should be repeated 
possibly also with other off-shell extensions. Moreover, 
other production mechanisms here omitted should be 
possibly included in the calculation, at least those mecha- 
nisms which proved to have relevant interference effects in 
NN collisions. But to implement the production mecha- 
nisms in a three-nucleon process is not a trivial task and 
needs to be done gradually. 

At present stage, where we believe that the importance 
of the off-shell effects in s-wave pion production from 
NN collisions has been fairly well established by various 
groups, it is clearly of great relevance to consider the 
consequences of such effects on more complex reactions. 
Herein we provide the results obtained when calculating 
off-shell effects in pd collisions. 



II. THEORY 

The production mechanisms are constructed starting 
from the following effective pion-nuclcon couplings: 
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(1) 
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The first term represents the gradient coupling to the 
isotopic axial current, while the second denotes the cou- 
pling to the isovector nucleonic current, and the last is 
the pion- nucleon coupling in the isoscalar channel. 

As is well known a good guess for the couphng 
constants can be obtained from chiral symmetry and 
PCAC, which constrain the three constants to be of the 
order 



AnXj/ml ^ 1/(4/2), 



ATrXQ/iriT^ ~ 0, 



(2) 



(3) 



(4) 



where gA (— 1.255) is the axial nucleonic normalization, 
and is the pion decay constant (~ 93.2 MeV). The first 
condition follows directly from the Goldberger-Trieman 
relation, while the last two are implied by the Weinberg- 
Tomozawa ones. Current phenomenological values for 
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fJ_MN/{4:n) can reach values as low as 0.0735 JlSt, .0749 
13 , or 0.076 |§, but also values around 0.081 have 
been considered acceptable. Some years ago common 
values were centered around 0.078-0.079 p2| , p3t . Simi- 
larly, from the pion-nucleon scattering lengths, is de- 
termined within the range 0.055-0.045, while the weaker 
isoscalar coupling has typically larger indetermination, 
ranging from 0.007 to -0.0013 [|,|. The isovector and 
isoscalar couplings, when combined with the axial irNN 
vertex, are the basic ingredients for the two-nucleon s- 
wave rescattering mechanisms, while the axial-current 
coupling alone forms the basis for the one-body produc- 
tion process. 

The matrix elements for the rescattering process re- 
quire an off-mass-shell extrapolation of the two constants 
Xj and Xq , since the rescattered pion can travel off-mass- 
shell. For Xq we consider the off-shell structure previ- 
ously employed in the pp — )■ pp?:" process by Hernandez 
and Oset (Ref. Q) which is based on a parametrization 
due to Hamilton M, namely 



(5) 



-(1 + e) + 



with TO^ = 550 MeV, a„ = 0.22 m^i, a^r = -0.233 m-\ 
and e = m-^jM. In the threshold limit, (g — p)^ ~ 
(qo — "t-tt)^ — q^, where q is the transferred 4- momentum 
between the two active nucleons. According to previous 
treatments of the IN 7r-production amplitude, the time- 
component qo is fixed around qo — while the space 
component q represents a loop variable and is integrated 
over. 

This form leads to an on-shell value of the order of 
0.007. The on-shell value derives from a cancellation 
between the short-range term, asr, and the tr-exchange 
term, am while off shell the cancellation occurs only par- 
tially and thus leads to the off-shell enhancement. The 
use of a fictions sigma-exchange model should not be 
considered a crucial aspect of the model, since similar 
forms (possibly summed over a "distribution" of masses 
rria) can be easily obtained also in theoretical approaches 
based on subtracted dispersion relations. The approaches 
based on subtracted dispersion relations such as TriV 
model of Ref. lead indeed to similar off-shell enhance- 
ment. 

For Xj the extrapolation can be accomplished by in- 
voking p-meson dominance and the related Riazuddin- 
Fayyazuddin-Kawarabayashi-Suzuki identity, which im- 
plies (on shell, for — > m^r) 
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with the corresponding off-shell extrapolation 
XjiqJ)^Xr9i{q,p) 



= X°r 



A2 

p 



{q-pf 



(7) 



where again we use (g — p)^ ~ {qo — rnT^Y ^ ^ i'^i the 
threshold limit. 

The production matrix-elements in the nonrelativistic 
2>N space with such couplings are the following: 



{iN\A^^\?,N,TT) = 
X ^(p' - p 



NN 



6-t-2e 



■cr2p[X2ji 



,(2 + ,)P^)^(q'-q+3P'') (8) 



for the one-body process. 



{3N\A^J^\3N,TT 



2i 



/■uNN'^'^Xq 



- xr- ^z^ v^NNiq)go{q,P) 



m2 + q2 - ' ■ 3 
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and 



{3N\Aj^\3N,7T) = V2*^^^^^4^ 



,,.^v^NN{q)gi{q,P) ,~ , ^ 



qo 



x^(q'-q+ip,) 



(10) 



for the two-body isoscalar and isovector rescattering, 
respectively. VT^NNiq) is the hadronic form factor of 
the ttNN vertex, whose structure is governed by the 
monopole-type cutoff A^. The momenta p and q are Ja- 
cobi momenta for nucleon 2 in the (2-1-3) center-of-mass 
(cm.), and nucleon 1 in the (1-1-2-1-3) cm., respectively, 
while is the pion momentum in the total cm. Sim- 
ilarly, p' and q' are the Jacobi momenta for the three 
nucleons in the no-pion case. Other relevant pion mo- 
menta are 
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and 



' 3(1 + e) 
(6 + 2e) 



6(2 + e) 



P. p' . 



(11) 
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In the actual calculation ranging from threshold up to the 
A resonance the on-shell couplings are further corrected 
by means of an Heitler-type (or iiT-matrix) form (77 is the 
pion momentum in pion-mass units 



A. 



1 Aq - 2A^ 



3 1 + 2iT]{XQ + Xj) 3l + 2i7j{Xo-2Xj) 



1 Aq + A, 

3 1 -t- 2ir;(Ao -f A^) 3 1 + 2ir]{XQ - 2Xj) 



1 Aq - 2A^ 



(13) 



(14) 
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This reduces the rescattering contributions at higher 
energies but the correction is uninfLuential in the thresh- 
old limit. On top of these processes, we have included 
also the two-body mechanism mediated by A rescatter- 
ing, 

{3N\Af\3N,7r) = ZlIlEA 



V'iVA(p',PA)52p[T2]^ 



Tcm + M-Ma+ Pa/Sa^a + q' I'lVA 
1, 



x5(q'-q+-P^), 



(15) 



since its contribution becomes soon important as the en- 
ergy increases. The intermediate A momentum is defined 
as 

In Eq. ( ]l5| ) /iA is the reduced mass of the A-iV sys- 
tem, while i^A is the reduced mass for the N-{AN) par- 
tition. Tc.m. is the cm. kinetic energy of the three 
nucleons in the initial state. The AiV transition po- 
tential is determined by the tt plus p exchange dia- 
grams, where the pseudoscalar meson provides the typ- 
ical longitudinal structure to the AA^ transition, i.e., 
(0-3 • q) (S'2 • q) (ts • T2 ) , while the vector meson generates 
the transversal contribution (0-3 x q)(>S'2 x q)(T3 • Tj). 
S2 (T2) is the spin (isospin) operator for the A ^ 
transition. For complete details on the employed transi- 
tion potential, and for other aspects connected with the 
isobar mechanism, such as, e.g., the detailed expression 
for the complex A mass A^a, we refer to a set of studies 
performed around the resonance region p^-p6[|. All am- 
plitudes, Eqs. (P),(p[),(p^), and ( p^ ) must be mu ltiplied 
in addition by the common factor l/v/(27r)3 2w^. More- 
over, taking into account Pauli identity, the full one-body 
mechanism results by multiplying Eq. (H) by the multi- 
plicity factor \/3(l + P), while the remaining two-body 
mechanisms are multiplied by 2-\/3(l + P), where P is 
the ternary permutator which commutes the 3A^ coordi- 
nates cyclically/anticyclically. Combining the P opera- 
tor with the given mechanisms in 3N partial waves is not 
a trivial task, and numerical treatment of the resulting 
amplitudes has been a challenge. 

The matrix elements are calculated between in and 
out nuclear states, where the out-state is specified by 
the three-nucleon bound-state wave function (plus a free 
pion wave) , and the incoming state is determined by the 
deuteron-nucleon asymptotic channel. For the 3A^ bound 
state we have taken the triton wave function as has been 
developed, tested and calculated in Ref. |2^. As two- 
body input for the three-nucleon equations we used a sep- 
arable representation of the Paris interaction. This 
form represents an analytic version of the PEST inter- 
action, originally constructed and applied by the Graz 
group pi. 



We have in addition calculated the relevance of the 
three-nucleon dynamics in the initial state (ISI) by solv- 
ing the Faddeev type Alt-Grassberger-Sandhas (AGS) 
equations . The AGS equations for neutron-deuteron 
scattering go over into effective two-body Lippmann- 



Schwinger equations |30| when representing the input 
two-body T-operators in separable form. The T is rep- 
resented in separable form using again the above men- 
tioned EST method. Applying the same technique to 
the TT absorption process, an integral equation of rather 
similar structure is obtained for the corresponding am- 
plitude. The only difference is that the driving term of 
the n-d scattering equation ( i.e., the particle-exchange 
diagram, the so-called "Z" diagram) is replaced by the 
off-shell extension of the plane- wave pion-absorption am- 
plitude. More details can be found in and refer- 
ences therein. 



III. RESULTS 



To exhibit the relevance of isoscalar off-shell effects 



for pd 



~t we have calculated the integral cross sec- 



tion near threshold up to the A region. The calcu- 
lated results are compared with a collection of data from 
Refs. ^ , ^J33| and others as explained in Fig. |l|. Prac- 
tically all the experimental data near threshold refer in 
fact to 77° production frompd collisions, and the compari- 
son has been made assuming isospin invariance and hence 
implying a factor of 2 between the two cross sections. In 
so doing we have avoided the need to include the effects 
of Coulomb distortions in the exit channel. Given the 
complexity of the reaction dynamics which depends upon 
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FIG. 1. Production cross section for the pd -K^t versus 
77. The dotted line contains the sole p-wave mechanisms. The 
dashed line includes also ttN s-wave rescattering mediated by 
p-exchange. The solid line considers in addition the isoscalar 
off-shell effects. The data are from Refs. [9,25,26]. The re- 
maining data ("others") have been extracted from a world 
collection as explained in Ref. [8]. 
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FIG. 2. The deuteron tensor analyzing power T20 for 
r; = 0.25. The lines show the same calculations as in Fig. 
1. The points are extracted from Ref. [9]. 
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FIG. 3. Effect of 3A'' initial-state correlation at = 0.25. 
Differential cross section {T20) on the top (bottom) panel. In 
both cases the solid line includes the SA^ ISI effects, while the 
dashed line has been obtained in plane-wave approximation. 
The points are extracted from Ref. [9]. 

several contributions, the isoscalar effects have been cal- 
culated on top of the other mechanisms we had consid- 
ered. As explained previously, the model includes also p- 



wave A rescattering, the one-body p-wave term, and the 
isovector p-exchange mechanisms. The relevant param- 
eters employed herein (cutoffs, coupling constants) have 
been tested previously against the pp —f Tr~^d process in 
Ref. p^. For the p-exchange model we use Xj — 0.045 
in Eq. iM). The results indicate that the modifications 
introduced by the isoscalar contributions are significant 
over the whole range considered, and that the effect is 
one of the most pronounced at threshold. 

Further evidences come from the results exhibited in 
Fig. ^ where the deuteron tensor analyzing power T2Q 
is shown. Details on the formalism for the calculation of 
polarization observables can be found in Ref. |34|. The 
points represent a best-fit to experimental data as given 
in Ref. |jl^. The trend of the data is reproduced once 
both the isovector and isoscalar terms are taken into ac- 
count. 

It is clearly important to assess the role of the ini- 
tial state correlations, since the three-body dynamics be- 
tween the nucleon and the deuteron could modify the 
whole picture and undermine the conclusions of this 
work. For this reason, we have calculated the ISI effects 
by solving the AGS equations for the 3N system using 
as input a separable representation of the Paris interac- 
tion. The same Faddeev-like technique herein employed 
has been applied previously to pion production at the A 
resonance in Ref. [|l5| . In Fig. |^ one can examine the role 
of the 3-/V dynamics in the initial state from the angular 
dependence of the unpolarizcd production cross section 
and from T20, for 77 = 0.25. The modifications intro- 
duced by the 3N dynamics are sizable, but the overall 
picture does not change drastically. In addition, the 3iV 
effects improve the angular dependence of both observ- 
ables, thus possibly confirming our findings about the 
importance of the isoscalar off-shell effects. 
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